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ABSTRACT 

Using mass-selected galaxy samples from deep multiwavelength data we investigate the incidence of close 
galaxy pairs between z = 0.4-2. Many such close pairs will eventually merge, and the pair fraction is therefore 
related to the merger rate. Over this redshift range the mean pair fraction is essentially constant (evolving 
as /p a i r ~ (1 + z)~° °' 6 ) with about 6 ± 1% of massive galaxies having a 1 : 4 or greater companion within 
30/i _1 kpc. Assuming the timescale over which pairs merge is not a strong function of redshift, this implies a 
similarly constant merger rate (per unit time) out to z = 2. Since about three times as much cosmic time passes 
at z < 1 as between z= 1-2, this implies that correspondingly more mergers occur in the low -redshift universe. 
When minor companions (1:10 mass ratio or greater) are included, the pair fraction increases to ~ 20% and 
still does not evolve strongly with redshift. We also use a rest-frame color criterion to select pairs containing 
only quiescent galaxies (major "dry merger" progenitors), and find them to be similarly rare and constant with 
4 - 7% of massive quiescent galaxies exhibiting a nearby companion. Thus, even though other studies find 
major mergers to be relatively uncommon since z = 1, our results suggest that few additional mergers occur in 
the 1 < z < 2 range and other mechanisms may be required to explain the mass and size growth of galaxies 
over this epoch. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — infrared: 
galaxies 



1. INTRODUCTION 

Most stellar mass in the universe was formed during 
a relatively brief chapter in its history, with the cosmic 
star-formation rate peaking around z ~ 2 (|Lillv et all [l 996; 
iMadau et al.|[T996t IHopkins & Beacoml l2006. and references 
therein). Cosmological models and simulations hold that local 
massive galaxies were hierarchically assembled, with mas- 
sive z = galaxies built from successive mergers of lower- 
mass constituents. Therefore, although many galaxies con- 
tinue to add mass through star formation at lower redshifts, as- 
sembly processes (such as major mergers and satellite accre- 
tion) should play an increasingly important role at late times 
(iGuo & White! 120091) . But recent observations reveal that 
the true picture is somewhat more nuanced: massive (M* > 
10 n M Q ) galaxies, many of which have ceased significant 
star formation, exist at leas t up to z ~ 2.5 dKriek et al.|[2 006, 
2008bt ICassata et alj|2008t [Stutz e t al J 120081: I Williams et all 
2009t iBrammer et al.ll2.009L 1201 lh . Although massive galax- 
ies should continue growing through mergers, a substantial 
fraction of their mass was in place at z = 2. 

Such hierarchical processes serve not only to build up the 
stellar masses of galaxies, but also affect their structural prop- 
erties (e.g. effective radius, surface density, and mass pro- 
file) and morphologies. It now appears that most quiescent 
galaxies are ext remely compact and den se for their mass 
at z ~ 2.3 (e.g. Ivan Dokkum et all 120081) . but comparably 
dense objects are p ractically nonexistent in the local universe 
dTavlor et al.ll2010b . Moreover, the average size of quiescent 
galaxies increases (and surface density wit hin the effective ra- 
dius decreases) smoothly with time (e.g. iFranx et alJ 120081 : 
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iDamianov etafl 120091: IWilliams et all 120 Id) . Gas-poor "dry 
mergers" have been proposed as a mechanism behind this 
growth, but the necessary evolution in the mass-size plane 
(much faster in size than mass) is difficult to achieve in equal- 
mass merging scenarios dBoylan-Kolchin et al] 12006). Sev- 
eral processes working in concert, inclu ding major mergers 
(e.g. lBell et aTJl2006l: IShankar et alj|201 lh . the quenching of 
progressiv ely more extended star -forming galaxies at lower 
redshifts (Ivan der Wei et al.ll2009l) . and the accretion of stellar 
mass in the outskirts of these compact "cores" (due t o mi- 
nor mergers; Bezans on et alJl2009l : IHopkins et alJl2009l) may 
therefore drive this size evolution. 

Mergers might thus play an important role in the buildup 
of massive galaxies, but observational determinations of their 
importance are still uncertain. At high redshifts mergers can 
be difficult to detect, thou gh clues may nonethe less appear in 
galaxies' morphologies (Conselice et al. 2003) and resolved 
velocity distributions dForster Schreiber et alJl2009h . Another 
promising technique is the use of close galaxy pairs as a proxy 
for mergers since these are easily detectable to high redshift 
and, on average, are expected to merge within a relatively 
short timescale (which can be calibrated using simulations; 
iBundy et al.ll2009hlLin et al.l201 U lkitzbichler & Whit d2008l 
KW08). 

Here we employ a deep 0.7 deg 2 multiwavelength survey 
to determine pair fractions between z = 0.4-2.0. First we 
present a brief summary of the data used in this sample, then 
the method for deriving galaxy pair fractions and the conver- 
sion to merger rates, and finally we discuss the implications 
of the pair fraction evolution for galaxy formation. AB mag- 
nitudes and cosmological parameters h = 0.7, Q,m = 0.3, and 
fl\ = 0.7 are used throughout. 

2. DATA AND DERIVED QUANTITIES 

For this study we use a catalog (Williams et al. 2012, in 
preparation) compiled from the UKIDSS Ultra-Deep Survey 
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(UPS: Lawrence et al.1 120071: IWarren et alj 120071) Data Re- 
lease 8 and supplementary data. The datasets and techniques 
employed for generating this catalog are sim i lar to those we 
used with the UDS DR1 in Willi ams et al.1 (120091): a brief 
summary follows. Source Extractor v2.5.0 dBertin & Arnoutsl 
1996) was run in dual-image mode to detect sources in 
the UDS DR8 K-band mosaic, measuring fluxes from a se- 
ries of deep optical/NIR images: u' from archival CFHT 
data, BVRi'z' from th e Subaru-XMM Deep Survey (SXDS; 
iFurusawa et all 120081) . and JHK from the UDS DR8, all 
convolved to the same point-spread function. Because of 
Spitzer's much larger point-response fu nction, we followed 
the deblending t echn ique described by lLabbe et"aT1 d2006l) 
and lWuvts et alj d2007l) to extract matched 3.6/xm and 4.5/um 
fluxes from deep IRAC imaging in the UDS field (SpUDS; PI 
J. Dunlop). Objects falling near bad pixels in the optical or 
near-IR images were excluded, as were those with no optical 
coverage, resulting in an effective image area of ~ 0.7deg 2 . 

From this updated catalog we calcula te photometric red- 
shifts with EAZY (Brammer et al. 2008 ), inte rpolate rest- 

120091) . and derive 
2009). The pho- 



frame colors with InterRest (Tavlor et al 
galaxy masses with FAST dKriek et a!T 
tometric redshifts agree to ~ 2% with spectroscopic red- 
shifts in the field; most catastrophic outliers are efficiently 
removed via a \ 2 cut. In calculating the masses we as- 
sumed a ISalpeterl (119551) initial mass funct ion (IMF), so- 
lar metallicity, and iBruzual & Chariot! (120031) stellar popula- 
tion models, and applied a fact or of -0 . 2 dex to the masses 
to bring them in line with a iRroupa d2001l) IMF (as in 
iFranx et all 120081: iWilliams et al.M201ot) . Although uncer- 
tainties on photometrically-derived quantities can be sub- 
stantial, in IR-selected galaxy samples stellar masses and 
mass-to-light ratios are relatively robust dKriek et alj [2008a; 
iMuzzin et al.l 120091) . The catalog has a 9 5% point-source 
compl eteness limit of K < 24.5; following Marchesi ni et al.l 
(2009), the corresponding mass completeness is estimated by 
scaling each galaxy to the mass it would have at a brighter 
K = 24.0 limit (to account for the fact that galaxies are not 
point sources); the 75th percentile of scaled galaxy masses 
thus provides an estimate of the 75% mass completeness limit, 
which we estimate as logM* = 9.8 at z ~ 2 for red galaxies 
(U-V KSt > 1.5). 

3. PAIR FRACTION MEASUREMENTS 

3.1. Major "wet" and "dry" merger candidates 

We searched the catalog for pairs of galaxies with projected 
transverse separations < 30/r 1 proper kpc, mass ratios of 1 : 4 
or greater, and require that the photometric redshifts of the 
galaxies in each pair lie within |zi -zi\ /(l +Zi) < 0.2 (where 
the subscripts 1 and 2 refer to the more and less massive 
galaxy in each pair, hereafter denoted "primary" and "sec- 
ondary" respectively). This large redshift separation ensures 
that few physically-assoc iated pairs are missed (a ccording to 
the analysis described in Ouadri & Williams 2010). From in- 
serting simulated pairs of point sources into the detection im- 
age, we find that SExtractor successfully deblends pairs with 
separations > l."2; we thus exclude those with < 13 kpc sep- 
arations (equivalent to 1 "5 at z ~ 2, to account for extended 
galaxy profiles). Only pairs containing one primary galaxy 
above Mi > 3.2 x 10 10 M Q and a secondary within the 1 : 4 
mass ratio constraint, i.e. 0.25 < M2/M1 < 1, are counted. 
The constraints on Mi and M2/M1 ensure that secondary 
galaxies are always above our \0gM2 > 9.8 completeness 



limit for red galaxies at z ~ 2. Finally, the constituent galax- 
ies within pairs are classified as either star-fo rming or quies- 
cent b ased on the rest-frame color criteria of Williams et al. 
(2009); each pair thus represents a candidate "dry merger" 
(i.e. quiescent-quiescent galaxy pair) or "wet merger" (con- 
taining at least one star-forming galaxy). 

Even with the relatively small angular separations and the 
requirement that galaxies within a pair be close in redshift, 
some apparent pairs are due to chance alignments. We correct 
for this by randomizing the positions of all galaxies over the 
survey area (while retaining the galaxies' masses, redshifts, 
and quiescent/star-forming classifications), and re-running 
the pair-finding algorithm on this rando mized catalog with 
the sa me selection parameters (see also lOuadri & Williams 
2010). This "mock catalog" process is repeated 20 times 
to reduce the uncertainty in the contamination rate, and the 
mean number of random pairs is subtracted from the total pair 
counts for each primary and secondary galaxy type. 

These corrected pair fraction measurements are listed in 
Table Q] and shown in Figure Q] We define the pair fraction 
fxy = Nxy/Nx as the fraction of "primary" galaxies of type X 
which have a less-massive "secondary" companion of type y; 
e.g., /sq is the fraction of log(M*/M Q ) > 10.5 star-forming 
galaxies with quiescent galaxy companions (whose masses 
are > 0.25M pr j mary ), counting each pair only once even if both 
members are above the primary mass threshold. The total 
pair fraction regardless of galaxy type, / a ii, a ii, is also given. 
Assuming these pairs merge at some point in the future, this 
fraction is related to the number of merger descendants - i.e. 
how many mergers a massive z = galaxy is likely to have 
undergone - as a function of redshift. 

3.2. Mass dependence and minor pairs 

While the pair sample in Table Q] makes full use of the 
galaxy catalog's dynamic range, it is also instructive to 
consider higher-mass primary galaxies: both to determine 
whether the pair frac t ion de pends on mass (e.g. Bund v et"aT] 
2009; Bernardi et al.l l20Tlb . and to investigate lower-mass 
companions which may represent minor mergers. The proce- 
dure described above is therefore repeated for primary galax- 
ies with log(Mi/M ) > 10.8 and companions within 1:4 and 
1:10 mass ratios; these fractions are shown in Figured Since 
star-forming galaxies are rare at this mass, only pairs with 
quiescent primary galaxies are included in this plot. Over 
all redshifts the 1:4 pair fractions are similar to those at 
log(M i /M Q ) > 10.5 shown in FigureHJ this is not surprising, 
since these samples are not widely separated in mass and the 
uncertainties are large. Unfortunately, the dynamic range of 
the catalog does not allow a detailed analysis of mass depen- 
dence due to the relative rarity of very massive galaxies and 
the need to maintain a "secondary" sample above the survey 
completeness limit. 

When minor companions are included, the pair fraction in- 
creases substantially: while only a few percent of these mas- 
sive galaxies have a close neighboring galaxy of comparable 
mass, 15-20% have companions within a 1 : 10 ratio (Figure|2l 
bottom panel). Since high-mass galaxies are more likely to be 
quiescent, this increased pair fraction is largely due to star- 
forming companions; however, with the inclusion of lower- 
mass secondary galaxies there even appears to be a marginal 
increase in the "dry pair" fraction at each redshift. Although 
the "minor pair fraction" is substantially enhanced, this does 
not necessarily suggest a high minor merger rate since satel- 
lite mass and dynamical friction timescale are expected to be 
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inversely correlated (KW08). 

3.3. Inferred merger timescales and rates 

Galaxy pairs only represent merger candidates, as some ap- 
parent pairs will not coalesce even over many Gyr while oth- 
ers rapidly merge, depending on orbital properties, projection 
effects, and other internal factors. Thus, converting from pair 
fractions to merger rates requires an average timescale, 

t merg = ^Vpair /N merger ( 1 ) 

where is the number of close pairs within a given redshift 
bin and N me i ger is the number of mergers per unit time within 
the same galaxy population. One theoretical parametrization 
of this timescale was derived from the Millennium Simula- 
tion by KW08 as a function of galaxy mass, redshift, and 
maximum projected separation. Specifically, their formula for 
pairs with projected separations Ar < 30/r'kpc, M2/M1 > 
1/4, and photometric redshifts gives a mean merger timescale 
of 

(?mer g (/r 1 Myr))" I/2 = 0.0189-9.47 x 10 _4 z 

+6.71 x 10 -3 [log(M*//r 1 M )-lO] (2) 

In the KW08 formalism the same mass limit is imposed 
on both primary and secondary galaxies, while our sam- 
ple includes lower-mass companions down to the survey 
limit. The resulting KW08 pair mass ratios depend on pri- 
mary galaxy mass (i.e., closer to equal-mass near the survey 
limit). For consistency, we thus repeat the above exercise 
with their selection method. Two mass thresholds are applied: 
log(M*/M Q ) > 10.5 as before, and log(M*/M ) > 10.0 to 
take full advantage of the survey. These fractions, shown 
in the top panels of Figure [3] are somewhat lower than in 
Figure Q] where secondary galaxies below the primary mass 
threshold are included. 

Figure [3] (bottom) shows the specific merger rates, calcu- 
lated with eq. (2), from z = 0.4-2 at the two mass thresh- 
olds. Adopting the KW08 timescale, mergers are relatively 
rare: only about 0.5-1% of massive quiescent galaxies merge 
with quiescent companions each Gyr, and fewer than 10% of 
galaxies over this redshift range will undergo major mergers. 
This is due both to the rarity of massive galaxy pairs and to 
the long effective merger timescales from KW08, about 2-3.5 
Gyr depending on mass and redshift. Although the KW08 
formalism is most applicable to our specific pair selection pa- 
rameters, a variety of timescales hav e bee n employed in the 
literature; this is discussed further in $4.11 

4. DISCUSSION 

4.1. Connecting pairs to mergers 

Close pairs of galaxies are easily identified and detectable 
to high redshift, and are therefore in principle a robust way of 
identifying systems that may merge within a relatively short 
timescale. But it is less straightforward to determine this 
timescale and convert the measured pair fractions to merger 
rates. 

Several estimat es of the merger timescale have been used 
in previous work. IB ell et al. (2006) assume that galaxy pairs 
merge within roughly one orbital time, in their case ~ 0.4 Gyr, 
while the KW08 estimate is nearly an order of magnitude 
larger for the galaxy masses considered here. These funda- 
mentally chang e the interp r etation of the measured pair frac- 
tions: with the B ell et al.l (12006b timescale, major mergers 



play a significant role in the assembly of massive galaxies 
over time; assuming KW08, only about 20% of massive qui- 
escent galaxies have undergone major mergers since z = 2, 
with ~ 2/3 of this occurring at z < 1. This disagreement is 
largely a result of the KW08 analysis including physically- 
associated galaxy pairs which are at relatively large real dis- 
tances despite having close projected separations, and there- 
fore merge only after a long period (or not at all). Nonethe- 
less, if the merger timescale isn't a strong function of redshift, 
the unchanging pair fraction we measure reflects a similarly 
constant merger rate since z = 2. 

4.2. Comparison to previous work 

The pair fractions de rived here are in bro ad agreement with 
previous work at z < 1 . iBundy et al.l (120091) find a low (~ 4%), 
non-evolving fraction of massive galaxies in "major pairs" at 
z < 1 .2, with 10 1 1 Mq galaxies more likely to have close com- 
panions than those with 10 I0 M Q . Given their higher mass 
limit and smaller search radius (20/i _1 kpc), our z < 1 mea- 
surements appear to be in agreement: a total pair fraction 
of ~ 4-8%, depending on whether the primary galaxies are 
star-forming or quiescent, and no strong evolution in the frac- 
tion. In addition, we confirm their reported higher incidence 
of "dry pairs" at lower redshifts, simply due to the coincident 
increase in the number density of massiv e quiescent galaxies. 
Between 0.4 < z < 0.8. lBell etai] (12006) report a pair fraction 
of 5 ± 1 % for galaxies above 2.5 x 10 10 M Q , also in agreement 
with our pair fraction measurement at the same redshift. 

One commo n theme in these studies is the rarity of dry 
mergers: even iBell et all d2006l) . with their short assumed 
timescale, find that only ~ 50% of m assive galaxie s have 
undergone major mergers since z = 0.8; iBundv et all d2009l) 
estimate 30% at the high-mass (10 u Mq) end (but see 
Padilla et al. 201 1). In the ir analysis of the environmental de- 
pendence of merger rate, iLin et al.l (1201 ll) estimate a some- 
what higher rate in high-density environments: 1 .2 ± 0.3 ma- 
jor dry mergers per galaxy since z = 1, perhaps not surpris- 
ing since such en vironments harbor a larger fraction of quies - 
cent galaxies (e.g. Kauff mann et al.l2004tlOuadri et al.l201 ll) . 
However, these only account for 38 ± 10% of the mass ac- 
cretion of these galaxies. If the merging timescale doesn't 
vary strongly with redshift, our results suggest that the major 
merger rate at 1 < z < 2 is comparable to that at z < 1 ; since 
twice as much cosmic time passes in the latter epoch than the 
former, this in turn implies that most major mergers occur be- 
low z ~ 1 . 

4.3. Major vs. minor mergers 

Even when short timescales are assumed, dry mergers from 
z = 2 to the present occur perhaps once or twice per galaxy 
at most. Since strong size and mass growth a re nonethe- 
less seen over this same redshift interval (e.g. iFranx et al] 
120081: IWilliams et al.ll20lol) . it appears that major mergers are 
not the primary driver behind the observed evolution. In- 
deed, if major dry mergers were the primary driver behind the 
smooth evolution in galaxy sizes and surface densities seen 
over z = - 2, a much larger number would be required to 
eliminate the comp act quiescent galaxy population by z = 
dTaylor et al.ll20Ta> . Even if mergers are more common than 
the KW08 timescale implies, they still may not account for 
the evolution of the mass-size relation since major mergers 
increase galaxies' masses and sizes at similar rates. 

Another explanation is that galaxies undergo many minor 
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mergers or accrete low-mass satellites. iNaab et alj (2009) 
find that minor mergers are more efficient (per unit secondary 
galaxy mass) at increasing galaxy radii than major merg- 
ers. This accretion scenario is also attractive because it de- 
pends only weakly (if at all) on whether or not the cen- 
tral galaxy is forming stars, and might also partially explain 
the size growth obse rved in massive star-forming galaxies 
dWilliams et al.ll20lol) . As shown in Figure [2] the inclusion 
of minor companions increases the pair fraction to 15-20% 
(see also Lopez-Sani uan et al.ll201 ll) . However, the merging 
timescale of these minor galaxies is not well-understood: on 
one hand their dynamical friction timescales are longer than 
more massive satellites, but tidal effects may hasten the incor- 
poration of minor satellites into a massive galaxy's extended 
halo. 

5. CONCLUSIONS 

From deep mass-selected samples, we have performed the 
first analysis of major and minor galaxy pairs out to z = 2, 
distinguishing between "dry" and "wet" merger candidates. 
"Dry" pairs are relatively rare: only ~ 3 — 7% of massive 
quiescent galaxies have a quiescent companion of compara- 
ble mass (with an average companion mass 40% that of the 



primary). This fraction increases slightly at lower redshifts, 
but the total pair fraction (including both star-forming and 
quiescent galaxies) remains essentially constant, evolving as 
/aii.aii ~ (1 +z)~ 0,4±0,6 . Minor companions are significantly 
more common, increasing to 15-20% for > 1 : 10 mass ra- 
tios. However, significant uncertainties remain when convert- 
ing pair fractions to merger rates. Nonetheless, while galaxies 
are often assumed to assemble rapidly at high z, the constant 
pair fraction suggests that, in fact, massive galaxies undergo 
more mergers at z < 1 than z= 1-2. Unless the effective 
merger timescale is much shorter than the values assumed 
here, such mergers represent a rare and stochastic process that 
is unlikely to make more than a marginal contribution to the 
smooth mass and size growth seen since z = 2. 
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TABLE 1 

Pair fractions in the UKIDSS-UDS 



-'■phol 


Nq 


N s 


/Oq 


/Qs 


fsq 


fss 


/all,all 


0.4<z<0.8 


270 


77 


0.066 ±0.017 


0.013 ±0.009 


0.035 ±0.026 


-0.001 ±0.013 


0.069 ±0.016 


0.8 <z < 1.2 


495 


163 


0.046 ±0.011 


0.016 ±0.008 


0.046 ±0.019 


0.015 ±0.014 


0.062 ±0.012 
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NOTE. — Uppercase subscripts denote the "primary" galaxy type: quiescent (Q) or star-forming (S); lowercase denote the type of the 
secondary galaxy within a pair. Hence "/q s " is the fraction of massive quiescent galaxies with a star- forming companion. Here all primary 
galaxies have masses above logMi > 10.5, secondaries are above logM2 > 9.9, and a mass ratio constraint of 0.25 < M2/M1 < 1 has been 
imposed. 
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FIG. 1. — Evolution of the massive galaxy pair fraction, defined as the fraction of massive "primary" galaxies (log(Mi/Mo) > 10.5) which have close 
companions within a 1:4 mass ratio (i.e. 0.25 < MijMy < 1; the average mass ratio is Mi/M\ ~ 0.4). The top and bottom panels show this fraction for 
quiescent and star-forming primary galaxies respectively. Red, blue, and grey points denote the fractions of quiescent, star-forming, and total secondary galaxies, 
respectively; thus, the red points in the top panel show the prevalence of major "dry merger" progenitors. The total pair fraction is remarkably constant over this 
redshift range. 
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FIG. 2. — Comparison of pair fractions with two different mass ratios (M7/M1 > 1/4 and Mt/M[ > 1/ 10) as a function of redshift. Here all "primary" galaxies 
are quiescent and 0.3 dex higher in mass than in Figure[T] black points show the fraction of primary galaxies with any close companions within the given mass 
ratio, while red points denote quiescent secondary galaxies ("dry mergers") only. In both mass ranges the pair fraction is constant or weakly declining with 
redshift. 
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FIG. 3. — Top panels: Pair fractions of quiescent primary galaxies determined using the same mass threshold for both primary and secondary galaxies. This is 
the selection method used in the theoretical merger timescale calibration developed by KW08. Bottom panels: Merger rates computed from the measured pair 
fractions and KW08 timescales. In all panels, red points denote quiescent-quiescent galaxy pairs ("dry mergers") while black points include both star-forming 
and quiescent secondary galaxies. The left and right panels show the same calculation for mass limits of logA/* /Mq = 10.0 and 10.5 respectively. Note the 
fractions in the right panel differ from Figure [T] due to the higher mass limit imposed on the secondary galaxies; in both cases the inferred merger rates are 
essentially constant and quite low. 



